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Abstract 

Quadratic curvature corrections to Einstein-Hilbert action lead in gen- 
eral to higher-order equations of motion, which can induced instability of 
some unperturbed solutions of General Relativity. We study conditions 
for stability of de Sitter cosmological solution. We argue that simple form 
of this condition known for FRW background in 3 -I- 1 dimensions changes 
seriously if at least one of these two assumptions is violated. In the present 
paper the stability conditions for de Sitter solution have been found for 
multidimensional FRW background and for Bianchi I metrics in 3 + 1 
dimensions. 

1 Introduction 

Recently theories of higher-order gravity have become a matter of intensive in- 
vestigations. Some of them are motivated by string theory fT], other proposals 
have been put forward in order to explain present accelerated expansion of the 
Universe (see [21 [3] and references therein) . One of the main features of these 
theories is the fact that they generally lead to equations of motion which con- 
tain higher-order derivatives. In 3 -I- 1 dimensions all higher-curvature terms in 
Lagrangian which keep the equations of motion to be the same order differential 
equations as in General Relativity vanish p[j . In string gravity the second order 
correction contains dilaton scalar field in addition to the gravity sector. These 
corrections has a particular structure (a product of Gauss - Bonnet term and 
a function of a dilaton field) which prevents the resulting equations of motion 
from being of higher-order equation. However, these equations become higher- 
derivative starting from the 3-d order corrections (see, for example, [5|.) 

As the number of time derivatives in the equations of motions increases in 
comparison with the situation in classical General relativity, the dimensional- 
ity of the phase space of the corresponding dynamical system increases as well. 
This may result in a dynamical behavior qualitatively different from that known 



in GR: some important GR solutions can loose their stability (this effect from 
a general point of view have been described in [6]). Though higher-derivative 
terms may be small on classical solution, perturbations from higher-order cor- 
rections can grow in direction of "additional" dimensions of the phase space, 
and a point which is stable in GR may become unstable in higher-order the- 
ory. This situation has already been described in 4-th order string gravity for 
Schwarzschild and de Sitter solutions [8J, as well as in cosmological scenar- 
ios with vacuum polarization [9]. The other proposal of this type is so called 
modified gravity with corrections in the form of function of the curvature f{R) 
[TOl [TT] which is also not free from instabilities [HI [H] (see, however, pH] and 
recent rewiews |15[ 116)). 

Recent understanding that modern string theories may admit many (ac- 
tually, very many) solutions in which 4-D part is the (metastable) de Sitter 
space-time with a positive effective cosmological constant ([iTl [l^ and many 
following papers) renew interest in the investigation of the stability of de Sit- 
ter solutions in gravity theories containing higher order powers of the Riemann 
tensor in their effective Lagrangian. These terms were long known to arise from 
quantum-gravitational corrections to the Einstein gravity (either the vacuum 
one, or interacting with matter quantum fields). Presumably, some of them 
may remain in the string theory in spite of many cancellations due to its numer- 
ous symmetries. Many particular forms of second order curvature terms have 
been studied in connection to possible instabilities induced by them |19 [ [2 m [2T|. 
This work have been also extended to i?" gravity \22\ I23j. In the present paper 
we study stability of de Sitter solution with respect to second order curvature 
corrections in their general form. 

2 FRW Universe in 3 + 1 dimension 

A general form of second order curvature corrections is 

S2^ J V^iaR^^'^^Rmm + ^3R'''R^k + iR^) d^x. (1) 
However, in 3 -I- 1 dimensions the Gauss-Bonnet term 

/-^ D riiklrn o a j^ik d i d2 

LiO — H rCikim — 4-n iXik + it 

does not contribute to equations of motion, and, then, we can rewrite (1) in the 
form 

^2 = y 7^(07?*"" i?,fei™ + iiR'^R.k + 7i?' - olGB) d^x. (2) 

Introducing B — /3-f 4q!, C — 7 — a we get the following action for the theory 
with a cosmological constant term 

5" / V^iR + BR'^R,k + CR^ - A) d-^x. (3) 
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For the Friedmann-Robertson-Walker metrics 



y,k = diag{-n{tf,a{t)^,a{tf, a{tf). 



(4) 



we have 



p2 
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Varying the action (3) with respect to n and setting n = 1 we obtain the 
following analog of Friedmann equation (we denote D — 12B + 36C): 

6H'^-A + 2DH^{H + H^)-3DH^-D{H + H^f + 2DH{H + SHH + H^) = 0, 

(5) 

where H = is the Hubble parameter. 

Note that the Friedmann equations, being algebraic in GR, becomes a differ- 
ential one in the theory with second order curvature corrections when D ^ 0. In 
the D = case the Friedmann equation does not change at all, and we exclude 
this case from the further analysis. 

For stability studies this equation should be written in the form of a system 
of two first order equations: 



H = F, 

F = -3HF 



2DH 



(A - QH^ + 2DH'^ + DF^) = f. 



(6) 



The form of de Sitter solution {H = 0,F = 0) give us the equations for 
corresponding stable points: 



1 



2DH, 



■iA-6Hi + 2DHfl)=0, 







its solution is Hq = iy ^. It is remarkable, that this de Sitter solution is 
exactly the same as in the pure GR, so higher-order corrections in the theory 
under investigation do not shift the location of the fixed point. However, its 
stability can be affected by the higher-order terms. 
Linearizing (jH we obtain 



H = F, 

F = (m] F- 



where 



OH 



-3H 



F 



(7) 



-3F - 3ff2 



5/ 

dF ■ 
A - 6H^ + 2DH^ + DF^ 



2DH^ 



12H + 8DH^ 
2DH ■ 
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Two eigenvalues of this system are 



For stability of the solution it is necessary that both eigenvalues have nega- 
tive real parts. Substituting the de Sitter solution we get that in a stable point 



As ( ) < in an expanding Universe, negativity of eigenvalues requires 



(^) ^ other hand, we have from (7) that in the de Sitter point 

df \ _ -6 

which means that the condition for stability of the de Sitter solution is £> = 
12(a + /3 + 37) > 0. 

In is possible to consider this result from a different point of view. Instead 
of Riemann tensor, we can express second order corrections through Weil tensor 

2 , CAr c^ d2 , M 13 Ar , 1 Ar2 \ d2 



where N is the dimensionality of space-time and the curvature scalar: 

5" = J V^{R + BC"'''"'C\kim+CR^-A)d^x, (8) 

where B = 2B = 2a+ ^(3 and 

C=l{B + 3C) = ^{a + p + 37) = £»/36. 

The Weil tensor vanishes on the Friedmann metrics, and the general quadratic 
corrections reduce to the i?^-tcrm with the known result: the de Sitter solution 
is stable if the coefficient before i?^-term in (8) is positive. 



3 Multidimensional FRW Universe 

It is instructive to compare this result with the situation in higher dimensional 
theory, where the Gauss-Bonnet term does contribute to dynamical equations. 
To illustrate importance of this property, we consider the theory with second 
order curvature corrections, which we write down in the form 

S = j V^{R + A&^^^Cikim + BE^^Rik + CR"" - A) rf^ar, (9) 

As usual, Cikim = on the Friedmann metrics, and we have two independent 
parameters of the theory B and C. The explicit forms of second order corrections 
on the Friedmann metrics are 
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R = 



[N-l) 



2-~2 + (iV-2) — 



(10) 



[N-l) 



"2 .... . n . 2 ■■ • 2 • • S -41 

N—-2N +N— — +2{N-2) — ^-2 -+{N-2f^ 



i?2 



'iklm 



2{N~1) Ld'^ 



naa , , a'* 

2 4— + 2 — -+ TV- 2 - 
a"^ naa n"^ a^ a* 



(11) 
(12) 



Varying ^ with respect to n, setting n = 1 and solving the equation of 
motion with respect to the second derivative we get 

H = 2mN-i)[B'N\4CiN-i)\ [iN - IKN 2)H^ - A + B{{N lf{N A)H^ 
-N{N - 1)H^ + 2N{N - i fHH^} + C{N{N ~ 1)'^{N - 4)i?4 - 4(iV - 1)^H^ 



+8{N- = 

(13) 

It should be noted that if BN + 4C{N — 1) = the resulting equation of 
motion does not contain H and H terms, which means that a new de Sitter 
solution is trivially stable. This case should be excluded from our analysis. 

The de Sitter fixed points can be found from the equation 



{N - 1)^{N - 4)(B + CN)H^ + (iV - l)(iV - 2)H^ - A = 0, 
and have the form 



-{N - 2) ± y^{N - 2)2 + 4A(iV - 4)(B + CN) 
2{N - 1){N - 4){B + CN) 



(14) 



(15) 



Comparing these results with those obtained in Sec. 2 we can see two impor- 
tant differences from the previously studied case of 3 + 1 dimensions: 

• The de Sitter points are shifted in comparison with classical point = 
K/{{N - 1){N - 2)). It happens due to the first term in eq.([l4l) which 
vanishes in iV = 4 dimensions. 

• There are two different de Sitter solutions, corresponding to two positive 
roots of l|15p . Each of them has its own zone of existence, so depending 
on B and C we can have two or one de sitter solutions, or the situation 
when de Sitter solutions are absent. 

The point with (+) sign in (fT5| has a regular behavior for B ^ 0, C ^ (it 
tends to unperturbed de Sitter solution), the behavior of the point with (— ) 
sign is singular. The latter point can not be an analog of any de Sitter point, 
existing in the pure GR, and exists even in the case of A = 0. 
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For stability analysis we should rewrite ifTSj) in the form of dynamical system 
(7). Using the same notations as in the previous section, we get 



( 



5/ 
dF 



) 







{N^l)Ho, 



and, for the same reason, the condition for de Sitter stability is equivalent to 
negativity of (^) • The eq. (fT3l) gives 



(H) =-i- 2H(A.-i)iB^+4C(^-i)l [2(^ - - 2)F + B{A{N lf{N - A)H^ 
+'iN{N - ifHH} + C{AN{N - lf{N - 'i)H^ + 16{N - l)'^HH}]. 



We get that stability conditions depend on the sign of the combination E = 
BN + 4:C{N —1). If it is positive, the regular de Sitter point (if exists) is always 
stable, while the singular de Sitter point is always unstable, independently of A. 
If it is negative, the regular point is unstable, and the singular point is stable. 

The conditions for existence of these de Sitter points can be easily derived 
from (15). Two conditions which should be satisfied are: {N — 2)^ + 4A(A^ — 
4)(B + CN) > and > 0. Corresponding zones on {B,C) plane are shown 
in Fig.l. 

Combining these results with the stability analysis we can distinguish the 
following zones in the parameter space {B,C) (see Fig. 2): 

• for negative E: 1) Zone with one unstable regular de Sitter (NOR), 2) zone 
with one stable (singular) and one unstable (regular) de Sitter (ROSK), 
3) zone without de Sitter solutions(if S'M). 

• for positive E: 1) Zone with one stable regular de Sitter (PON), 2) zone 
with one stable (regular) and one unstable (singular) de Sitter (POSL), 
3) zone without de Sitter solutions {LSM). 



4 Bianchi I Universe in 3 + 1 dimensions 

The situation become less simple even in 3 + 1 dimensions if we allow anisotropic 
perturbations of the de Sitter solution, which means that the Weil tensor has a 
non-zero contribution to the equations of motion 0. 

Consider a flat homogeneous anisotropic Universe with the Bianchi I metrics 
in the form 

^Main results of this section have been found independently using another set of variables 
in [25]. 




In the de Sitter point we get 




df \ _^VW 



2)2 + 4A(A^-4)(B + CiV) 
BN + 4:C{N- 1) 



(16) 
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Figure 1: Zones of de Sitter solutions for the action (9). A singular de Sitter 
exists in double shaded band zone between line a [B+CN > —{N—2Y/ {AK{N— 
4))] and b [i? + C7V > 0], a regular de Sitter exists in both double shaded band 
and shaded zone above the line b. 



g,fc = diag{-n{tf, a{tf ,b(t? , c{tf )- (17) 

It should be noted that in the presence of general quadratic curvature cor- 
rections the diagonal form of the metrics (fTZl) is not the general one [24J. We 
restrict ourself by this special form for simpHcity. We show that even in this less 
general case we have additional restrictions for stability conditions of de Sitter 
solution. 
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Figure 2: Zones of de Sitter stability for the action (9). The Hnes LK and PR 
are the same as in Fig.l; the Hne MN is the E = line. The regular de Sitter 
(if exists) is stable above this line, the singular de Sitter is stable below this line 

Starting from the action (3) we have three equations of motion 

2Hb + 2He + 2Hl + 2Hl + 2HMb -L + 2B(i/f ' + Sij^ + 2Hi^'> + + H^Ht 
+3HcH, + ZHl + Hf'^ + -Ht + hH.HaHa + HaH.Hb - HlHaH, - H^HaHb 
+HbH^Hc + AHcHbHb + + mcH^ + AHcHa + 2HcHb + HaHb + 2HbHc 
-HlHa + HlHb - HlH, - HlHb + HlH, - H^Ha + 2H^H^ + HaA - 2HlH^ 
+AH^Ha + 'iHlHb + 2,HbHb + '^HbH^ + iHaHb + HaHl + iHaHc + HaHl 
+H^Hc + 3HbHc + 2HbH^ + 211^11^ — HaHbHb + 4:HbHcHc — HaHcHc 
+HbHaHc + 5HbHaHa + AHMbHa) + 4C(2i7i'' - HlHl + QHMc + 
+bHl -H^ + + 2iff ^ + ml + 2H^^^ + ml - 2HbHlH^ + 2HlHb - 2H^Hb 
-2HlH^ + 2HlH^ - H^Hl + mlHl + 2HaHe + 2HaHb + ^^H^b + AHbHa 
+QHlHb + QHbHb - 2HlHa + 4HaHa + AHcHa + GHcH^ + AHaHb + 2HaHl 
+6HbHc + AHbHl + AHlH^ + AHbHc + AHaH, + 2i?„,ij,i?, + 2HaHbHb 

+2HbHaHa + AHcHbHa + 2HaHl + SHcHbHb + 2HcHaHa + 2HbHaHc 

+8HcHbHc + 2HcHaHb) = 0, 

(18) 
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where Ha = a/a, Hh = b/b, He = c/c, two others can be obtained by transmu- 
tation a ^ c and a ^ b. 

Substituting the de Sitter solution {Ha = Hf, = He = Hq =Const), we 
get that the de Sitter solution has its unperturbed form Hq = i/A/6. After 
evaluating the corresponding system of equation of motion, we find the following 
eigenvalues 



Mi,2,3 = -3//o, (19) 

_ -I8CH0 - 6BH0 ± 2 v'81C2i72 + MBCH^ + 9B^H^ -&C-2B 
^^'^ ~ 2(6C + 2B) ' 

(20) 

-WHo ± y/33B^H'^ +4B + 9QBCH^ . , 

M6,7,8,9 = ^ • (21) 

In expanding Universe we have always /ii,2,3 < 0. The condition /i4,5 < 
leads to the already known result B + 3C > 0. However, we have now additional 
restrictions arising from the condition /Lt6,7,8,9 < 0. It leads to AiJ + 4CA+ 1 < 
if S > and A\B\ - AC A - 1 < if B < 0. Note that these stability conditions 
now depend not only on coefficients before the higher-order terms in the action 
(3), but also on the value of the cosmological constant A. These additional 
restrictions arise from instability with respect to anisotropic perturbations of the 
de Sitter metrics. It is interesting that independently of A the whole quadrant 
B > 0, C > 0, which always satisfies the stabfiity conditions in a FRW Universe, 
is excluded in anisotropic case (see Fig. 3) 

In the representation (8) we have the following conditions 4CA < BA/6 — 1 
for S > and 4C'A > \B\A/6 - 1 for B < 0. (see Fig.4). 



5 Conclusions 

We have studied the influence of second order curvature corrections in their gen- 
eral form on stability of de Sitter solution. Our results for the 3-1-1 dimensional 
space-time and for higher-dimensional cases differs significantly. The reason of 
this difference is non-zero contribution from the Gauss-Bonnet combination for 
higher-dimensional space-times. In 3 -|- 1 dimension the de Sitter solution al- 
ways exists, though it may be stable or unstable depending on particular form 
of the second order corrections. The stability condition takes a very simple and 
independent of the value of cosmological constant A form for isotropic Universe, 
though become more complicated in anisotropic Universe. 

In multidimensional Universe even in isotropic case (the only considered 
in the present paper) there are several possibilities. Depending on the action 
wc may have two de Sitter solutions ( in this case one solution is stable and 
one solution is unstable), one de Sitter solution (which can be either stable or 
unstable), and there are theories with no de Sitter solutions at all. 

There results can be useful in constructions of cosmological scenarios taking 
into account higher-derivative terms. Instability of a regular de Sitter solutions 
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1 



Figure 3: Zones of de Sitter stability for the action (3) (shaded) The de Sitter 
solutions in zone below the line 1 are unstable on isotropic background, the 
solutions in the unstable zone above the lines 1 and 2 are unstable due to 
anisotropic perturbations. 

can damage the standard inflationary scenario when close to de Sitter infla- 
tionary stage is produced by an effective cosmological constant originated in a 
matter sector of the theory. This opens a principal possibility to rule out some 
theories of this type from the observational point of view. On the other hand, 
unstable singular de Sitter solution can be used in alternative approach to infla- 
tion |26j. Our results show that quadratic curvature corrections in Lagrangian 
can produce such kind of solutions only in multidimensional Universe. 
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Figure 4: Zones of de Sitter stability for the action (8) (shaded) The unstable 
zone in the upper half-plane appears due to anisotropic perturbations. 
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